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Abstract. In this short note we consider the empirical distribution of the eigenvalues of the sum of 
d independent Haar distributed Unitary matrices, and show that this sequence of measures converge 
to the Brown measure of free sum of d unitary operators. As a consequence of our proof we improve 
the uniform bounded condition on the imaginary part of the Stieltjes transform in [4, Theorem 1]. 

1. Introduction 

Random matrix models can be broadly divided into two regimes — Hermitian matrix models, 
and non- Hermitian matrix models. Many methods, like the method of moments, and the Stieltjes 
transform work well in the hermitian regime, and provide a good understanding of the behavior 
of the spectrum for many different matrix models. However both of these methods fail in the 
non-hermitian regime, and understanding the behavior of the spectral measure in any given model 
is much more difficult. 

In an article by Girko [3], a nice scheme was proposed to find the limiting spectral distribution of 
non-hermitian matrix model. However carrying out this scheme, even in the simplest model, where 
all the entries are i.i.d. with mean and finite variance, is extremely challenging. After a long 
sequence of partial results (see references in [2]), the circular law conjecture, for the i.i.d. case, was 
recently established by Tao and Vu [10] in full generality. Barring this simple model very few results 
are known in the non-hermitian regime. For example, nothing is known about the spectrum of 
random oriented d- regular graphs. It was recently conjectured in [2] that the empirical distribution 
of the adjacency matrix of a random oriented ci-regular graph sampled uniformly converges to a 
measure (i on the complex plane, which has a density 

1 d\d-l) 

Vr(d2- \ z \2)2\\z\<Vd} 

with respect to the lebesgue measure on C. This conjecture is due to the observation that fi is the 
Brown measure of the free sum of d unitary operators (see [6, Example 5.5]), and motivated us to 
consider a simpler problem, about the sum of d independent unitary matrices. For the latter, we 
prove the following: 

Theorem 1.1. For any positive integer d, the empirical distribution of the eigenvalues of the sum 
of d independent, Haar distributed unitary matrices, converges weakly, in probability, to the Brown 
measure of free sum of d unitary operators. 
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It is worthy to note that if {U^}f =1 is a sequence of independent Haar distributed unitary 
matrices, then it converges in ^-moments (see [9] for a definition) to {ui}f =1 , where Uj's are free 
Haar unitary operators (see [1, Theorem 5.4.10]). However Brown measure being discontinuous, 
convergence in ★-moments does not necessarily imply the convergence of Brown measures 1 (see [9, 
§2.6]). Nevertheless, it is proved in [9, Theorem 6] that if the original matrices are perturbed by 
adding small Gaussian (of unknown variance), then the Brown measures do converge. But removing 
the Gaussian, and even identifying the correct variance, is often a very hard task. For example 
in [5, Proposition 7, Corollary 8] it is shown that for a particular ensemble adding polynomially 
vanishing Gaussian is not sufficient to regulraize the Brown measure. Also generalizing the notion of 
convergence in ^-moments to the notion of convergence in distribution of traffics [7] , does not help 
to remove the Gaussian. We circumvent these obstacles, and in Theorem 1.1 we essentialy show 
that the original matrices are smooth enough, and therefore there is no need to add any Gaussian 
in this case to get convergence of the corresponding Brown measure. 

To prove Theorem 1.1 we adopt the techniques developed in [4]. There one relies on the uniform 
boundedness of the imaginary part of the Stieltjes transform (see [4, Theorem 1, Eqn. 3]). Here, 
the expected empirical distribution of the eigenvalues of \U n — zl n \ has unbounded density (see 
Lemma 3.1), and therefore the imaginary part of its Stieltjes transform is unbounded. Nevertheless, 
in this short note we show how to control the unbounded regions so that we achieve the desired 
result, namely the convergence of the integral of the logarithm near zero for Lebesgue almost every 
z. It is not hard to realize that, using similar argument, same result as in Theorem 1.1 can be 
obtained, if we replace unitary matrices by Haar distributed orthogonal matrices. 
In [4, Theorem 1] it is shown that under assumptions [4, Eqn. 1-3] the empirical distribution of 
the eigenvalues of U n T n converges, is a single ring, and is related to the Brown measure of the 
corresponding limiting operator. In a recent paper by Rudelson, and Vershynin [8] it is shown that 
the condition [4, Eqn. 3] can always be removed (see [8, Corollary 1.3]). It can also be noted that 
[4, Eqn. 2] is a very rigid condition, and it prevents the existence of atoms in the support of the 
limiting distribution. For example, it can not include the simplest case where T n = I n (see [4, 
Remark 2]). But as a corollary of our proof, we improve the uniform boundedness condition of the 
Stieltjes transform (compare (1.2) with [4, Eqn. 2]) to obtain the following which can accommodate 
atoms: 

Corollary 1.2. Let empirical distribution of the eigenvalues of {T n } converges weakly to a proba- 
bility measure 0, compactly supported on K+. Assume further 

(1) There exists a constant M > so that 

lim P(||T n || > M) = 0. (1.1) 

n— >oo 

(2) There exists a set of isolated points {xi}f =l such that, for every e > 0, there exists < 
k,M < oo for which the following hold: 

{z : Qf(z) > n" K , \%(G Tn (z))\ > M} C {z : G uf =1 ( Xi -e, Xi + e)}. (1.2) 
If Q is not a Dirac measure, then the following hold. 

(a) The empirical distribution of the eigenvalues of A n := U n T n converges, in probability, to limiting 
probability measure 

(b) The measure jjla possesses a radially- symmetric density with respect to the Lebesgue measure 
on C, satisfying pa(z) = ^A z (f log \x\du z (x)), where A z denotes the Laplacian with respect 
to the variable z, v z := EH A| 2 |, A r = \[b r + 5- r ), and 6 is the symmetrized version of 0. 



For a matrix its Brown measure is its empirical distribution of the eigenvalues (see [9, Proposition 1]) 
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(c) The support of [Xa is single ring: There exists constants < a < b < oo so that 

suppfiA = {re 10 : a < r < b}. 

Further, a = if and only if f x~ 2 d@(x) = oo. 

This extension above can now include the cases where can have atoms, or unbounded density, 
as long as, at the finite n-level (1.2) satisfied. For example, if T n is diagonal, with ai proportion of 
them are Xi, for i = 1,2, ... ,k, then G is X{ with probability pi, and at the finite n-level (1.2) is 
satisfied. The case when T n = al n , for some a > 0, is immediate from Theorem 1.1. 

2. Proof of Theorem 1.1 

Theorem 1.1 is proved via Girko's method, and for completeness, we start with describing the 
key steps in Girko's method (for more details see [4]). To this end, given a matrix A n , let La„ 
denote the empirical distribution of its eigenvalues. For every z £ C, define 

(A n - zl n 

_{A n - Zl n )* 

and let be the empirical distribution of the eigenvalues of H z . Girko's method consists of: 

Step 1: Show that for (Lebesgue almost) every z E C, the measures converge weakly, in 
probability, to measure v z as n — > oo. 

Step 2: Justify that f log \x\v^{dx) — > J log \x\v z (dx). 

In the non-hermitian regime, the main (technical) challenge is to verify Step 2. Beyond that 
completing the proof of Theorem 1.1 requires few relatively standard additional steps (outlined in 
sequel). 

Turning to prove Step 1 and Step 2, for any z E C, and {U n }i<i<d i.i.d. unitary matrices, distributed 
according to Haar measure, define, 

l,z ._ \ {U n ~ Zl n 

(U n ~ Zl n )* 



H' 



V, 



l.z 



u. 



and 



rd,Z 



-1,2 



[U. 



d\* 



yd-l,z , 

r TV. 



"0 u d ~ 




0" 


^ n 




+ 


{U^* 



for d > 2. 



For any two compactly supported measures n± and /i2 on M, we denote /Uiffl/X2 to be free convolution 
(see [1, §5.3.3] for more details) of /Ui and fi2- 

Step 1 and Step 2 of Girko's method are taken care of by the following two lemmas respectively 
(whose proofs are deferred to Section 3). 

Lemma 2.1. Let 1,z is the symmetrized version the measure having density given by (3.1) with 
r = \z\, and for d>2, Q d ' z := 6 d_1 ' z EE Ai, where Ai = \{6~-i + h)- Then for d E N, the empirical 
measure L v d,z converges weakly, in probability, to Q d,z . 

Lemma 2.2. For Lebesgue almost every z E C, for every d > 2, 



J log \x\dL v d,z(x) 



\og\x\d® d > z {x), 



(2.1) 



in probability. Furthermore, for every d E N there exists {rj}|Lj (depending on d), such that for 
every e > and R > 0, whenever <f> is any smooth function supported on {z : \z\ < R}\ {z : 
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\z\ G (rj — e, r, + e)}, we have 



z) I log\x\dL v d,z(x)dm(z) 



log \x\dQ d,z (x)dm(z) 



(2.2) 



in probability. 



Proof of Theorem 1.1: First let us assume that (2.2) holds almost surely. Combining Lemma 2.2 
with [4, Eqn. 5], we get that for every R > and e > 0, and any smooth function ip, supported on 

(2.3) 



T d R<£ := 5(0, R)\A a e :={z:\z\< R}\ U^'i {z : M G {n -e,n + e)}, 



we have 



1 

~2^ 



Aifj(z) I log\x\d<d d > z (x) dm{z) 



Hz) 



log\x\de d ' z (x) dm(z):= I %j>(z)h(z)dm{z) , (2.4) 



where the integration by parts in the last but one step follows from tp being compactly supported 
and smooth. Now by [6, Proposition 3.5] we note that Q d,z is the symmetrized version of the law 
of \ui + U2 + • • • + Ufi — zl\, and from the definition of Brown measure it follows that h(z)dm(z) is 
the Brown measure of u\ + U2 + • • • + (see [6]). We further note that for d > 2, it has a density 
(namely, h(z)) with respect to the Lebesgue measure on C. Therefore by a standard argument it 
follows that for any open set GcC, for any R > 0, and e > 0, 



liminfL f/1+ ... +c/d (Gn(r£ e )°) > f 



h(z)dm(z), 



a.s. and consequently 



liminf Lrn , , n d(G) > lim lim / 



h(z)dm(z) 



Gn(r% )•> 



> lim 



h(z)dm(z) — lim / h(z)dm(z) = / h(z)dm(z), (2.5) 
'C7nB(0,R) £ ^°J(Af)° Jg 

almost surely. Since the same statement remains true for a countable collection of open sets, the 
proof completes by a use of Portmanteau theorem, if (2.2) holds a.s.. 

To prove the above when (2.2) holds in probability, we use the fact that the convergence in probabil- 
ity implies that for every subsequence, there exists a further subsequence such that a.s. convergence 
holds. To this end, for any given subsequence,we can extract a further subsequence (by cantor diag- 
onal argument) such that (2.5) holds for a countable collection of opens sets G C C, and therefore 
again by a use of Portmanteau theorem, the proof is completed for d > 2. 

When d = 1, from the proof of Lemma 2.1 it is clear that the expected empirical measure of the 
eigenvalues of \U n — zl n \ is same as ' z , and therefore the expected Brown measure of U n is same 
as that of Haar unitary. Consequently the proof is completed by a use of a standard concentration 
inequality for Haar measures on unitary group (see [ I , Corollary 4.4.30]). □ 

3. Proofs of Lemma 2.1 and Lemma 2.2 

Since U^s are invariant under rotation, it is easy to see that Q d ' z is a function of \z\. Therefore 
the lemmas will be proved only for z = r > 0. Before we prove Lemma 2.1 we will need a preliminary 
result first: 
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Lemma 3.1. Let U n be Haar distributed. For any r > 0, let Lu nr be the expected empirical 
distribution of the eigenvalues of \U n — rl n \. Then Lu n r has a density 

2 x 

f r (x) = , |r-l|<x<r + l. (3.1) 

W vr^(x 2 -(r-l) 2 )((r + l) 2 -x 2 ) ' " V ' 

with respect to the Lebesgue measure. 

Proof. Note that to prove the above it is enough to show that the expected empirical distribution 
of the eigenvalues of (U n — rI n )(U n — rl n )* has the density 

gJx) = - — 1 =, (r-1) 2 <x< (r + 1) 2 . (3.2) 

vr x /( x -( r _l)2)(( r + i)2_ 3 .)' ^ > - -\ ) \ ) 

To this end first note that using rotational invariance of the Haar measure, we have Tr{[/^'}] = 
Tr{([/*) fc }] = for any positive integer k. Thus 

- Tr {(U n + U*) k } = ( J for k even and otherwise. (3.3) 

L Un+u * = 2 cos 6 = e ie + e~ id , where 6 ~ Unif (0, 2ir). (3.4) 
Now noting that 

(U n - rI n )(U n - rl n )* = (1 + r 2 )I n - r(U n + U*), 
the result follows from the change of variable formula. □ 
Proof of Lemma 2.1: Let L, r d,r be the expected empirical distribution of the eigenvalues of V r f' r . 

The fact that I v i,r = © 1,r follows from Lemma 3.1. We proceed to show that L^ r d,r =^ Q d,r , 
by induction on d > 2, following the proof of [4, Lemma 10], setting there p = 1, and replacing 
T n by Vn~ l)1 ', to conclude (relying on the Schwinger-Dyson equation (3.5)) that the corresponding 
Stieltjes transforms converge. The fact that L„d,r Q d,r , in probability, follows from the standard 
concentration inequality for Haar measures (see [1, Corollary 4.4.30]). □ 

Lemma 2.1 completes the proof of Step 1 of Girko's method. Turning to Step 2 for d > 1, and 
z E C + , define 

G d n r (z)=E\±Tr(zL n -V n d n- 



E 

Therefore 



2n 



and 

,r 1 



Gg» = E[— Tr [u d (zL n - F^)" 1 } 



Now imitating the steps leading from [4, Eqn. 30] to [4, Eqn. 38], we obtain the following finite 
n-level Schwinger-Dyson equation: 

G d /(z) = G d -^(^ d n (z)) - 0(n,z,ip d n (z)), (3.5) 

(analogous to [4, Eqn. 38]) where 

ri(z):=z- . (3.6) 

l + 2G d d :(z) 

(analogous to [4, Eqn. 37]), for every z\,Z2 £ C + 

Q(n,z 1 ,z 2 ) 



0(n,z 1 ,z 2 ) - 



l + 2G^( 2l ) 
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and 



0(n, zi,z 2 ) =E 



O 







' 1 1 




-E 








_2n . 










-E 




\2n 




I2n J 



d{ Zl i n - v n d n~\z 2 i n - vt 1 *)-^* 



i 



^ 2 |9(z 2 )|9(zi) 2 (9(zi) M)J- (3 ' 7) 

We note that the Schwinger-Dyson equation (3.5) is valid whenever the denominator in (3.6) does 
not vanish. In particular, it holds when $s(z) is large. However our aim is to extend the Schwinger- 
Dyson equation to a region very near to the real line. To do so we need the following fine estimates 
on the terms appearing in (3.5). To this end similarly as in [4, Eqn. 34] we further obtain, 

(3.8) 



(Griz))' = 2G#»(1 + 2G£») - Oi (*>*). 



d.r 



where 



Oi(n,z) 



4E 
O 

















-E 






-E 




\2n 




fin - 






.2n . 



d(zl n - V r 



d.r' 



-Ijjd 



Thus, 



n 2 %(z) 2 (%(z) A 1) 



G*» = \ ( - 1 + V 1 + + 40x(n, z) 



(3.9) 



d.r i 



f» w "4\~ 1T V 1T,lu " w "r^U"-)*^! (3.10) 
where the branch of the square root is uniquely determined by analyticity, and the behavior of G^f 

and Gn T when \z\ — >■ oo . Furthermore as in [4, Lemma 11, Lemma 12] we get: 

For every positive integer d>2, and < r < R, there exists constants C\, C2,C^ depending only 

on d and R such that for all z G C + , with Ss(z) > Cin -1 / 3 , and all large n, 

|1 + 2G*»| >C 2 [9(z) 3 Al], (3.11) 

and for all z G C + such that ^s(z) > C^n~ 1 ^, 

*Mz)) > 9(*)/2. (3.12) 

Thus it shows that the finite n- level Schwinger-Dyson equation (3.5) extends to the region {z : 
^s(z) > Can -1 / 4 }, for all large n. 

The proof of Lemma 2.2 consists of controlling the Stieltjes transform of G n ' r (z) inductively in d, 
where Lemma 3.2 provides the basis d = 1 of the induction, and Lemma 3.3 the inductive step. 

Lemma 3.2. For any e > and r > there exists some constant C , such that 

{z: \^G l /{z)\ >Ce~ 2 } C L := E + in G C + : n G (0,£ 2 ),£G ( ± (1 ± r) - 2s, ±(1 ± r) + 2ej J . 
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Proof. First note that for n > 0, 

\%G l /(E + i v )\ = J)^ f r ( x ) dx + / Jl^ f r ( x ) dx 

J\x-E\>y/fj \ x ~ +V J\x-E\<^rj \ x — J^) +V 



< 1 + 

< 1 +7T 



1 f V 

SUP / r (x) / =r r - = 



r/.r 



sup / r (x) . (3.13) 

x:\x- E\< y/rj 

Denoting T e to be the union of open intervals of radius e around the four points ±1 ± r; from (3.1) 
it follows that for any e > 0, 

sup{/ r (x)} = Cie" 1 , 

for some constant C\, independent of choice of r. Thus from (3.13) it follows that 

sup l^G 1 /' (E + < Ce~\ (3.14) 

Now noting that 

{(E,r)):Ee T c 2e , r) G (0, e 2 )} C {(E, rf) : (E — y/rj, E + yjrj) G T^}, 

and 

sup |$G£ r (£ + ir/)| < e" 2 , 

the proof follows from (3.14). □ 
In Lemma 3.2 we note that < 3Gn ,r (z) is not uniformly bounded, and therefore results similar to 
[4, Lemma 13] is not possible in our set-up. Instead in Lemma 3.3 we show that if we have some 
knowledge about the regions where Q(G n ' r (z)) blows up, using (3.5) we can identify the regions 
where ^s(G r ^' 1 ' r (z)) blows up, which is good enough to establish Lemma 2.2, and consequently 
Theorem 1.1. 

Lemma 3.3. Fix any d G N, and some r G (0, R\. Given any e > assume that there exists 
constants K, M , and K , depending possibly only on d and R such that 

T d := {z : > n~ K , |3(G# p (z))| > M} C {z : R(z) G Vf r }, (3.15) 

where rf' r is the union of e intervals around ±m ± r, for m = 1,2, . . . , K . 

Then there exists constants k\,M\, and K±, depending possibly only on d and R such that, 

T d+l := {z : > n~ K \\Q(G d n +1 ' r (z))\ > M x } C {z : G T^ 1 '*}, (3.16) 

where rf +1 ' r is the union of e intervals around ±m ± r, m = 1, 2, . . . , K\. 

Proof. Enlarging M in (3.15) as needed, we assume without loss of generality 6M _1 < e 2 . Fixing 
any < K\ < k, and Mi = 2M let z G 1^+1. Now note that by (3.7) and (3.11), we get 

\0{n,z,^ +1 {z))\ < -rn < M, 

n 2 |9(^n +1 (^))|9(2) 2 (9(z) 4 Al) " 

whenever ^s(z) > C'n~ K ' for some C' < oo, n' > 0, and using (3.12) in the right most inequal- 
ity. Thus from (3.5) and (3.12), we get |9(Gn' r \^ +1 (z)))\ > M, and %{^ +1 (z)) > n~ K , and 
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consequently tp^ (z) G T^. Now for any z G C + define 



R{G d +^{z)) :- 



G d n +1 > r (z) 



Using (3.8) we note that 



G 



d+i,n 



(z) l-R 2 (G d n +1 > r (z)) =R{G^{z)) + 



d+l,r, 



G 



d+l,r 



(z)Oi(n,z) 



(l + 2G^)) : 



(3.17) 



(3.18) 



Now using (3.10) we further obtain that, 

2Gt +1 > r (z) 



R(G d n +1 ' r (z)) 



l + ^l + 4(C^ 1,r (s)) 2 + 40i(n,z) 
1 G d+1 ' r (z)Jl + 4(Gt +1 ' r (z)) 2 + 40i(n, z) 1 



G d+1 ' r (z) 



2 (Gi +1 ' r (z))^ + 1 (n,z) 2(^( z ))2 + 0l(nj , ) - (3J9) 

Noting that, from (3.9), Oi(n, z) = 0(n _1 ) for ^s(z) > n" 1 / 3 the rightmost term in the RHS of 
(3.19) is bounded by M _1 whenever |3?(<4 +1,T "(,z))| > M. Also note that if both %(z) > n" 1 / 3 , and 
|9(Gn +1 ' r (z))| > M, then for any choice of the branch of the square root, 



1 G d+1,r (z)Jl + 4(G d n +1 > r (z))i + 40x(n,z) 



(G d n +1 ' r (z)y + 1 (n,z) 



< 2k 



|l + 4(G^ +1 ' r (z)) 2 + 40i(n,z) 



\A(G d n +1 ' r (z)y\ 



< 4. 



Therefore 

\R(G d ^(z))\ < 5, 

whenever ^(z) > ?i -1 / 3 , and |3 t (Gn +1 ' r (^))| > M. Furthermore, from (3.9) and (3.11), there exists 
constant k>2 > such that 

m > „-« =► gg^jgi(n^) = 0(n -i } 
Therefore, from (3.18) we get that for Q(z) > n~^ 2A1 / 3 \ and |9(G^ +1,r »)| > M, 



tf 2 (G^(z)) - 1 



This further shows that #(G^ +1,, ») G B(l,e) U J3(-l,e) 5(z,e) := {z' : |jk - z'| < e} . 

Now collecting all the arguments above we see that if we choose k\ small enough (in particular any 
K\ < K A K2 A 1/3 will do) , and M\ = 2M then for any z G r^ + i, we have ifj d+1 (z) G T^, and from 
(3.6) it follows that 



R(«) = K(^ +1 (z)) +5R( J R(G^ 1 ^(z))) G rg- 1 -*-. 



□ 



Thus the proof is completed. 
Now we are ready to prove Lemma 2.2. 

Proof of Lemma 2.2: First note that by Lemma 3.2, Lemma 3.3, and using induction, for every 
e > 0, and R > fixed, we get that for d G N, and r G (0, R], there exists constants k,M, and i'T, 
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depending only on d and R, such that for every 

z G {z : > n- K ,R(z) $ T d £ ' r }, 

\Q(Gi r (z))\ < M, where Tf r is the union of e intervals around ±m ± r, m = 0, 1, 2, . . . , K. Now 
we note that G Te' 1 " r G U^ 1 (rj — e,rj + e) for some sequence {rj}^ 1 . Therefore using [4, 
Lemma 15], and imitating the steps leading to [4, Eqn. 49], for any a G [1,2], and for any e' > 0, 
we obtain 



E 



r \ log x\ a dL v ^, z (x)] <Ce'\log{e')\ a , (3.20) 

J n~ K 

whenever z G (defined in (2.3)), and the constant C depends only on R. Now from [8, Theorem 
1.1] we know that there exists constants < c±, C2 < oo such that for every z G C, 

P(*mm(V^*) < t) < t Cl n C2 , 

and thus upon choosing a 5 > 02/ci, there exists a <5' such that for all large n, 



E 







logd a (2L. /d , 2 (x) <y. (3.21) 



We further note that without loss of generality we can assume k < 5 and, hence imitating the steps 
in the proof of [4, Proposition 14 (i)] we obtain (2.1). 

Now to prove (2.2) we first note that the bounds in (3.21), (3.20) are uniform for z G T^ £ , and 
then the proof is completed by adapting the steps in [ I, Proposition 14(h)] . □ 

Proof of Corollary 1.2: The main step of the proof is to show that for Lebesgue almost every z, 
logarithm is uniformly integrable with respect to the empirical distribution of the eigenvalues of 
\U n T n — zl n \. As shown in [4], this problem is equivalent to that of Y* := T n + p(U n + U*), where 
p = \z\. We note that for this latter problem, one key-step is to show that the imaginary part of the 
Stieltjes transform of the expected empirical distribution of the eigenvalues of T n + p{U n + £/*) is 
uniformly bounded on the imaginary axis, provided the imaginary part is bigger than n~ K , for some 
k > 0, for Lebesgue almost every p (see [4, Proposition 14 (i)]). In [4] this was shown under the 
assumption [4, Eqn. 2]. Here we show that same can be achieved under the modified assumption 
(1.2). 

To this end note that the Schwinger-Dyson equation obtained here, and that of [4], are similar, it 
is not hard to realize that, adapting the proof of Lemma 3.3, from (1.2), we obtain that for every 
e > 0, and R > there exists some constants M\, and k\ depending only on R and e such that, 
for every p G (0, R], 

{z : 9f(s) > n- Kl ,|9(G T „ +pM) (z))| > Mi} C {z : K(z) G uf =1 ( Xi ±p-2e,x i ±p + 2e)}. 

Now following the proof of Lemma 2.2, and [ J, Corollary 1.3] we obtain that for Lebesgue almost 
every z 

log \x\dv*(x) — > / log \x\du z (x), 



in probability, where is the empirical measure of the eigenvalues of Y* . Moreover the same proof 
shows that, for any small e > (any e < mm{xi : X{ 7^ 0} will do), any R > 0, and any smooth 
function 0, supported on B(0, R)\ Vjf =l {z : \z\ G (x{ — e, X { + e)}, we have 

J (f){z) j log \x\du^ l {x)dm{z) — > j 4>{z) j log \x\dv z (x)dm(z). 
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Since not a Dirac measure, from [6, Theorem 4.4], and [4, Remark 8] we note that \ia has a 
density with respect to the Lebesgue measure on C. Hence, combining the proof of Theorem 1.1, 
and [4, Theorem 1], we obtain the desired result. □ 

Acknowledgements. We thank Ofer Zeitouni for suggesting to look at the sum of d unitary 
matrices, and pointing our attention to preliminary versions of [2, 8]. 
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